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Abstract: An important line of research is the investigation of the laws 
of random variables known as Dirichlet means as discussed in Cifarelli and 
Regazzini ((7J) . However there is not much information on inter-relationships 
between different Dirichlet means. Here we introduce two distributional 
operations, which consist of multiplying a mean functional by an inde- 
pendent beta random variable and an operation involving an exponential 
change of measure. These operations identify relationships between differ- 
ent means and their densities. This allows one to use the often considerable 
analytic work to obtain results for one Dirichlet mean to obtain results for 
an entire family of otherwise seemingly unrelated Dirichlet means. Addi- 
tionally, it allows one to obtain explicit densities for the related class of 
random variables that have generalized gamma convolution distributions, 
and the finite-dimensional distribution of their associated Levy processes. 
This has implications in, for instance, the explicit description of Bayesian 
nonparametric prior and posterior models, and more generally in a variety 
of applications in probability and statistics involving Levy processes. We 
demonstrate how the technique applies to several interesting examples 
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1. Introduction 

In this work we present two distributional operations which identify relationships 
between seemingly different classes of random variables which are representable 
as linear functionals of a Dirichlet process, otherwise known as Dirichlet means. 
Specifically the first operation consists of multiplication of a Dirichlet mean 
by an independent beta random variable and the second operation involves an 
exponential change of measure to the density of a related infinitely divisible 
random variable having a generalized gamma convolution distribution (GGC). 
This latter operation is often referred to in the statistical literature as exponen- 
tial tilting or in mathematical finance as an Esscher transform. We believe our 
results add a significant component to the foundational work of Cifarelli and 
Regazzini (§0). In particular, our results allow one to use the often considerable 
analytic work to obtain results for one Dirichlet mean to obtain results for an 
entire family of otherwise seemingly unrelated mean functionals. It also allows 
one to obtain explicit densities for the related class of infinitely divisible random 
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variables which are generalized gamma convolutions, and the finite-dimensional 
distribution of their associated Levy processes, (see Bertoin (jl|) for the formali- 
ties of general Levy processes). The importance of this latter statement is that 
Levy processes now commonly appear in variety of applications in probability 
and statistics. A detailed summary and outline of our results may be found in 
section 11.21 Some background information, and notation, on Dirichlet proce- 
ses and Dirichlet means, their connection with GGC random variables, recent 
references and some motivation for our work is given in the next section. 



1.1. Background and motivation 

Let X be a non-negative random variable with cumulative distribution function 
Fx- Note furthermore for a measurable set C, we use the notation Fx(C) to 
mean the probability that X is in C. One may define a Dirichlet process random 
probability measure, see flih and (14; [HI), say Pg, on [0, oo) with total mass 



parameter 9 and prior parameter Fx, via its finite dimensional distribution as 
follows; for any disjoint partition on [0, oo), say (Ci, . . . ,Ck), the distribution 
of the random vector (Pg(Ci), . . . ,Pg(Ck)) is a fc-variate Dirichet distribution 
with parameters (0Fx(Ci), . . . , 9Fx{Ck))- Hence for each C, 

Pg{C) = / I(xeC)Pg(dx) 
JO 

has a beta distribution with parameters (9Fx{C), 9(1 — Fx{C))). Equivalently 
setting 9Fx(Ci) = 9i for i = 1, . . . , k, 

(Pe(C 1 ),...,P e (C k ))^(^;i = l > . 

where (GgJ are independent random variables with gamma(0i, 1) distributions 
and Gg — Gg 1 + • • • + Gg k has a gamma(#, 1) distribution. This means that 
one can define the Dirichlet process via the normalization of an independent 
increment gamma process on [0, oo), say 7e(-)j as 

Pe{-)- ^ 



7fl([0,oo)) 



where je(Ci) = Gg i and whose almost surely finite total random mass is 7e([0, oo)) 
Gg. A very important aspect of this construction is the fact that Gg is indepen- 
dent of p£, and hence any functional of Pg. This is a natural generalization of 
Lukacs'(|33f) characterization of beta and gamma random variables, whose work 
is fundamental to what is now referred to as the beta-gamma algebra, (for more 
on this, see Chaumont and Yor ((0), section 4.2)). See also Emery and Yor (11) 



for some interesting relationships between gamma processes, Dirichlet processes 
and Brownian bridges. 

These simple representations and other nice features of the Dirichlet process 



have, since the important work of Ferguson (|l4t Il5f ). contributed greatly to the 
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relevance and practical utility of the field of Bayesian non and semi-parametric 
statistics. Naturally, owing to the ubiquity of the gamma and beta random 
variables, the Dirichlet process also arises in other areas. One of the more inter- 
esting, and we believe quite important, topics related to the Dirichlet process is 
the study of the laws of random variables called Dirichlet mean functionals, or 
simply Dirichlet means, which we denote as 



/>oo 

M e (F x )= / xPe(dx), 
Jo 



initiated in the works of Cifarelli and Regazzini (@;0). In (0) the authors ob- 
tained an important identity for the Cauchy-Stieltjes transform of order 9. This 
identity is often referred to as the Markov-Krcin identity as can be seen in 
for example, Diaconis and Kemperman (fiol ). Kerov |26l ) and Vershik, Yor and 



Tsilevich (|37l). where these authors highlight its importance to, for instance, 
the study of the Markov moment problem, continued fraction theory and ex- 
ponential representation of analytic functions. This identity is later called the 
Cifarelli-Regazzini identity in (fig). Cifarelli and Regazzini (0), owing to their 
primary interest, used this identity to then obtain explicit density and cdf for- 
mulae for Mg(Fx)- The density formulae may be seen as Abel type transforms 
and hence do not always have simple forms, although we stress that they are 
still useful for some analytic calculations. The general exception is the case of 
6 = 1 which has a nice form. Some examples of works that have proceeded along 



these lines are Cifarelli and Melilli (|8j), Regazzini, Guglielmi and di Nunno (|35f ). 
Regazzini, Lijoi and Priinster|36h. Hjort and Ongaro (fl8h. Lijoi and Regazz- 
ini and Epifani, Guglielmi and Melilli (12; 13)). Moreover, the recent work 



of Bertoin, Fujita, Roynette and Yor (]2|) and James, Lijoi and Prunster (|23l ) (see 



also (2m) which is a preliminary version of this work) show that the study of 
mean functionals is relevant to the analysis of phenomena related to Bessel and 
Brownian processes. In fact the work of James, Lijoi and Priinster (f23h iden- 
tifies many new explicit examples of Dirichlet means which have interesting 
interpretations. 

Related to these last points, Lijoi and Regazzini (30) have highlighted a close 
connection to the theory of generalized gamma convolutions(see (|3|)). Specifi- 
cally, it is known that a rich sub-class of random variables having generalized 
gamma convolutions (GGC) distributions may be represented as 

/>oc 

T e = GeMg(Fx) = / x le (dx). (1.1) 



We call these random variables GGC(#, Fx)- In additon we see from (jl.ip that 
Tq is a random variable derived from a weighted gamma process, and hence 
the calculus discussed in Lo (f3lh and Lo and Weng (pjih applies. In general 
GGC random variables are an important class of infinitely divisible random 
variables whose properties have been extensively studied by (3) and others. We 
note further that although we have written a GGC(9,Fx) random variable as 
GgMg(Fx) this representation is not unique and in fact it is quite rare to sec 
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Tg represented in this way. We will show that one can in fact exploit this non- 
uniqueness to obtain explicit densities for Tg even when it is not so easy to do 
so for Mg(Fx)- While the representation GgMg(Fx) is not unique it helps one 
to understand the relationship between the Laplace transform of Tg and the 
Cauchy-Stieltjes transform of order 9 of Mg(Fx), which indeed characterizes 
respectively the law of Tg and Mg(Fx). Specifically, using the independence 
property of Gg and Mg(Fx), leads to, for A > 0, 

E[e- ATe ] = E[(l + \Mg(F x ))' e ] = e-<^ (A) (1.2) 

where 

roc 

i/j Fx (X) = / log(l + Xx)Fx(dx) = E[log(l + XX)}. (1.3) 
Jo 

is the Levy exponent of Tg. We note that Tg and Mg(Fx) exist if and only if 
^_f x (A) < oo for A > 0,(see for instance (0) and (0)). The expressions in (|1.2| 
equates with the the identity obtained by Cifarelli and Regazzini (0) , mentioned 
previously. 

Despite these interesting results, there is very little work on the relationship 
between different mean functionals. Suppose, for instance, that for each fixed 
value of 9 > 0, Mg(Fx) denotes a Dirichlet mean and (Mg(Fz );c > 0) denotes 
a collection of Dirichlet mean random variables indexed by a family of distribu- 
tions {Fz c ;c> 0). Then one can ask the question, for what choices of X and Z c 
are these mean functionals related, and in what sense? In particular, one may 
wish to know how their densities are related. The rationale here is that if such a 
relationship is established, then the effort that one puts forth to obtain results 
such as the explicit density of Mg(Fx), can be applied to an entire family of 
Dirichlet means (Mg(Fz c ); c > 0). Furthermore since Dirichlet means are asso- 
ciated with GGC random variables this would establish relationships between 
a GGC(9, Fx) random variable and a family of GGG(9, Fz c ) random variable 
Simple examples are of course the choices Z c = X + c and Z c = cX, which, due 
to the linearity properties of mean functionals, results easily in the identities in 
law 

Mg(F x+ c) = c + Mg(Fx) and Mg{F cX ) = cM e (F x ) 

Naturally, we are going to discuss more complex relationships, but with the 
same goal. That is, we will identify non-trivial relationships so that the often 
considerable efforts that one makes in the study of one mean functional Mg (Fx ) 
can be then used to obtain more easily results for other mean functionals, their 
corresponding GGC random variables and Levy processes. In this paper we will 
describe two such operations which we elaborate on in the next subsection. 

1.2. Outline and summary of results 

Section 11.31 reviews some of the existing formulae for the density and cdf of 
Dirichlet means. In Section [2l we will describe the operation of multiplying a 
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mean functional Mg a (Fx) by an independent beta random variable with pa- 
rameters (8(7,9(1 — a)), say, (3ea,e(i~a) where < a < 1. We call this op- 
eration beta scaling. Theorem 12.11 shows that the resulting random variable 
fieo,e(i-o)Me< T (Fx) is again a mean functional but now of order 9. In addi- 
tion, the GGC(#cr, Fx) random variable Gg r7 Mg r7 (Fx) is equivalently a GGC 
random variable of order 0. Now keeping in mind that tractable densities of 
mean functionals of order 9 = 1 are the easiest to obtain, Theorem 12.11 shows 
that by setting 9 = 1, the densities of the uncountable collection of random 
variables (/3 (Ti i_ cr M (T (Fx); < a < 1), are all mean functionals of order 6=1. 
Theorem 12.21 then shows that efforts used to calculate the explicit density of 
any one of these random variables, via the formulae of (0), lead to the explicit 
calculation of the densities of all of them. Additionally, Theorem 12 . 2 1 shows that 
the corresponding GGC random variables may all be expressed as GGC random 
variables of order 9 = 1, representable in distribution as G\{i a ^ a M a (Fx)- A 
key point here is that Theorem 12.21 gives a tractable density for ji a ^^ a M a (Fx) 
without requiring knowledge of the density of M a (Fx), which is usually ex- 
pressed in a complicated manner. These results also will yield some non-obvious 
integral identities. Furthermore, noting that a GGC(6*,Fx) random variable, 
Tg, is infinitely divisible, we associate it with an independent increment process 
(Ce(t)]t — 0) known as a subordinator ,(a non-decreasing non-negative Levy 
process), where for each fixed t, 

E [ e -ACe(t)] = E[e" ATet ] = e - te *pFxW. 

That is, marginally (e(l) = T$ and (g(t) = (et(l) = Tgf In addition, for 
s < t, (o(t) — (e(s) = (e(t — s) is independent of Co(s). We say that the pro- 
cess (C,e(t);t > 0) is a GGG(9,Fx) subordinator. Proposition 12.11 shows how 
Theorems 12.11 and 12. 2\ can be used to address the usually difficult problem 
of describing explicitly the densities of the finite-dimensional distribution of a 



subordinator (see (|27l )). This has implications in, for instance, the explicit de- 
scription of densities of Bayesian nonparametric prior and posterior models. 
But clearly is of wider interest in terms of the distribution theory of infinitely 
divisible random variables and associated processes. 

In Section^ we describe how the operation of exponentially titling the density 
of a GGC((9, Fx) random variable leads to a relationship between the densities of 
the mean functional Mg(Fx) and yet another family of mean functionals. This 
is summarized in Theorem 13.11 Section 13.11 then discusses a combination of the 
two operations. Proposition 13. II describes the density of beta scaled and tilted 
mean functionals of order 1. Using this, Proposition 13.21 describes a method to 
calculate a key quantity in the explicit description of the density and cdf of mean 
functionals. In section[4]we demonstrate how our results can be applied to extend 
and explain results related to two well known cases of Dirichlet mean functionals. 
However, we emphasize that Proposition 14. 1[ 14.21 and 14.31 are genuinely new 
results to the literature. More complex applications, which may be viewed as 
extensions of section I4T21 may be found in an unpublished preliminary version of 
this work in |2lh. We discuss and develop these further in James (22). Section 
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5 presents a more involved result relative to those in section 4, but which does 
not require a great deal of background material. Here we show how the results 
in section 2 are used to derive the finite dimensional distribution and related 
quantities of a class of subordinators recently studied in (0) . 

1 . 3. Preliminaries 

Suppose that X is a positive random variable with distribution Fx, and define 
the function 

/>oo 

$**(*) = / log(\t-x\)I(t^x)F x (dx)=E[log{\t-X\)l(t^X)]. 
Jo 

Furthermore, define 

Ae(t\F x ) = -sin(n6F x (t))e- e * F xW, 

TT 

where using a Lebesque-Stieltjes integral, F x {t) — f F x {dx). Cifarelli and 
Regazzini (|7j) (see also (|8|)), apply inversion formula to obtain the distributional 
formula for Mg(Fx) as follows. For all 9 > 0, the cdf can be expressed as 

(x-tf^AoitlFx^t (1.4) 

provided that 9Fx possesses no jumps of size greater than or equal to one. If we 
let £,8F X (') denote the density of Me(Fx), it takes its simplest form for 9=1, 
which is 

£ Fx (x) = A 1 (x\F x ) = -sm(nF x (x))e-^ x l (1.5) 
Density formulae for 9 > 1 are described as 

&f x (x) = (0-1) / {x~t) e - 2 A e {t\Fx)dt. (1.6) 
Jo 

An expression for the densit y, w hich holds for all 9 > 0, was recently obtained 
by James, Lijoi and Priinster (|23T ) as follows, 

Zsf x (x) = - [ {x-t) e - l d {t\F x )dt (1.7) 

where 



n 



deit\F x ) = ^sin{ir9Fx{t))e- e ^ . 
dt 

For additional formula, see iffi [351: l30h. 

Remark 1.1. Throughout for random variables R and X, when we write the 
product RX we will assume unless otherwise mentioned that R and X are in- 
dependent. This convention will also apply to the multiplication of the special 
random variables that are expressed as mean functionals. That is the product 
Mg(Fx)Mg(Fz) is understood to be a product of independent Dirichlet means. 
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Remark 1.2. Throughout we will be using the fact that if R is a gamma random 
variable, then the independent random variables R, X, Z satisfying RX = RZ 
imply that X = Z. This is true because gamma random variables are simplifiable. 
For precise meaning of this term and conditions, see Chaumont and Yor |j, sec. 
1.12 and 1.13). This fact also applies to the case where R is a positive stable 
random variable. 

2. Beta Scaling 

In this section we investigate the simple operation of multiplying a Dirichlet 
mean functional Me (Fx) by certain beta random variables. Note first that if 
M denotes an arbitrary positive random variable with density /m, then by 
elementary arguments it follows that the random variable W = f3 a j,M, where 
Pa,b is beta(a, b) independent of M, has density expressible as 

•^H = / fM(w/u)u a - 2 (l - uf-'du. 

r (o)r(6) J 

However it is only in special cases where the density fw can be expressed in even 
simpler terms. That is to say, it is not obvious how to carry out the integration. 
In the next results we show how remarkable simplifications can be achieved when 
M = Mg(Fx), in particular for the range < 6 < 1, and (3 a .b is a symmetric 
beta random variable. First we will need to introduce some additional notation. 
Let Y a denote a Bernoulli random variable with success probability < a < 1. 
Then if X is a random variable with distribution Fx, independent of Y a , it 
follows that the random variable XY a has distribution denoted as 

FxrAdx) = aF x (dx) + (1 - a)5 (dx), (2.1) 

and cdf 

FxyJx) = <tF x (x) + (1 - a)l(x > 0). (2.2) 
Hence, there exists the mean functional 

/•OO 

M 9 {F X Y a ) = \ yPe(dy) 
Jo 

where Pg(dy) denotes a Dirichlet process with parameters (9, Fxy a )- In addition 
we have for x > 0, 

§ FxVa (x) = Epogflz - XY a \)l{XY a ± x)\ = a$ Fx (x) + (1 - a) log(z). (2.3) 

When a — 1, Y a = 1 and hence Fxy x (x) — Fx(x). Let E a denote a set such that 
¥\Po(E a )\ — a. Now notice that every beta random variable, /3 a ,b, where a, b are 
arbitrary positive constants, can be represented as the simple mean functional, 

P e (E a ) ^ f3 eaM1 _ a) £ M e (F Ya ), 
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by choosing 

a 

a = and 9 = a + b. 

a + b 

We note however that there are other choices of Fx that will also yield beta 
random variables as mean functionals. Throughout we will use the convention 
that /3g t o :— 1, that is the case when a = 1. We now present our first result. 

Theorem 2.1. For 9 > and < a < 1, let ^e<r,e(i-<r) denote a beta ran- 
dom variable with parameters (9a, 0(1 — a)), independent of the mean functional 
M 9a (F x ). Then 

ft) Pe.,e{i-.)M 6 .{Fx) = M {F XYa ). 
(ii) Equivalents, M e {F Y „)M e<T {F x ) = M (F XY J. 

(ill) Gg a Mg a (F x ) = GgMe(F X Y a ) 

(iv) That is, GGC(9a, F x ) =GGC{6, F XY J. 

Proof. Since Mg(F Ycr ) = Pea,g(i-a) statements (i) and (ii) are equivalent. We 
proceed by first establishing (hi) and (iv). Note that using (|f .3[) . 

/>OC 

E[log(l + XXY a )] = (xE[log(l + XX)} =<j log(l + Xx)F x (dx). 

Jo 

Hence 

E r e -AG e M 9 (F X y <T )l _ c - e ° io S (l+\x)F x (dx) _ ^r-XG^M^Fx)] 

which means that GgMg(FxY CT ) = Gg a M 0a (Fx) 1 establishing statements (hi) 
and (iv). Now writing Gg a = Gg/Jeo-^i-o-)- It follows that 

GgMg(F X Y CT ) = Ggf3g a g( 1 _ a )M ga (Fx)- 

Hence f3g a ,e(i-a)Mea(Fx) = Mg{F X Y„), by the fact that gamma random vari- 
ables are simplihable. □ 

When 9 = 1, we obtain results for random variables j3 a .i^ a M„(Fx)- The 
symmetric beta random variables P a ,i-a arise in a variety of important contexts, 
and are often referred to as generalized arcsine laws with density expressible as 

Now using (HU) and (H]), let C(F X ) = {x : F x (x) > 0}, then for x > 0, 

M*F*yA*)) = { S[n ^f- Fx {^ l X l C r% X \ ( 2 - 4 ) 
v " K " \ sm(7r(l - a)), iixfC(F x )- 

Note also that sin(7r[l — Fx(x)}) = sin(7rFx(x)). The next result yields another 
surprising property of these random variables. 
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Theorem 2.2. Consider the setting in the Theorem \2.1i Then when 9 = 1, 
it follows that for each fixed < a < 1, the random variable M\(Fxy„) = 
/3 C r ! i_ (T M .(Fx) /las density 

£ Fxy<r (as) = — sin(7rF xr „ (x))e^*^ W /or z > 0, (2.5) 

TT 

specified by {2.1$. Since GGG(a,Fx) = GGG(l,FxY cr ), this implies that the ran- 



dom variable G a M a (Fx) = G\M\(FxYa) has density 

i r°° 

9a, Fx (x) = - e-iy^- 2 sin(TTFxYAy))e- a ^ iv) dy (2.6) 

TT Jo 

Proof. Since Mi(FxY a ) = l3 a ,i-aM a (Fx), the density is of the form (L5|) . for 
each hxed cr G (0, 1]. Furthermore we use the identity in (|2.3j) . □ 

Remark 2.1. It is worthwhile to mention that transforming to the random vari- 
able l//3 CT ,i_ CT , H2.5\) is equivalent to the otherwise not obvious integral identity, 

shuW) r~ ^ Fx ^v) dy = ^ sin{jTFxY (jB))e — (.). 

n Ji \V ~ i) ^ 

TTus /eads to interesting results when the density £ ct f x (#) ftas a known form. On 
the other hand, we see that one does not need the explicit density of M a (Fx ) to 

obtain the density of Mi(Fxy„) = 0<r i-a-M a (Fx)- In fact, owing to our goal of 
yielding simple densities for many Dirichlet means from one mean, we see that 
the effort to calculate the density of M\(Fxy^), for each < cr < 1, is no more 
than what is needed to calculate the density of M\ (Fx ) ■ 

We now see how this translates into the usually difficult problem of describing 
explicitly the density of the finite-dimensional distribution of a subordinator. In 
the next result we use the notation Co(C) to mean I (,s e C) (g(ds). 

Proposition 2.1. Let ((e(t):t < 1/9) denote a GGC(0, Fx) subordinator on 
[0,1/9]. Furthermore let (Gi,...,Cf.) denote an arbitrary disjoint partition of 
the interval [0, 1/9]. Then the finite- dimensional distribution (Ce(Ci), ■ ■ • , Ce(Cfc)) 
has a joint density 

k 

l[g^,F x (x l ), (2.7) 

i=l 

where each o~i = 9\Ci\ > and XL=i °~i = 1- The density g ai ,F x * s given by \2.b}) . 
That is, Co(Ci) = G\Mi(FxY a .) and are independent for i ~ l,...,k, where 
M 1 {F XYai ) = (S^x-^M^Fx) has density 

-x"*- 1 sin(irF X Y (x))e- a '^^ {x) . 

TT 
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Proof. First, since (C%, . . . ,Ck) partitions the interval [0,1/0], it follows that 
their sizes satisfy < \Ck\ < 1/6 and J2i=i | Cfe j = 1/6. Since is a sub- 
ordinator the independence of the Ce(Cj) is a consequence of its independent 
increment property. In fact these are essentially equivalent statements. Hence, 
we can isolate each Ce(Ci). It follows that for each i the Laplace transform is 
given by 

^[e-AC^c,)] = e ~e\c^ Fx (\) = e -<ni> Fx (\)^ 

which shows that each Ce(Ci) is GGC(<7i, Fx) for < <Xj < 1. Hence the result 
follows from Theorem 12.21 □ 



3. Exponential Tilting/Esscher Transform 



In this section we describe how the operation of exponential tilting of the density 
of a GGC(#, Fx) random variable leads to a non-trivial relationship between a 
mean functional determined by Fx and 6, and an entire family of mean func- 
tional indexed by an arbitrary constant c > 0. Additionally this will identify a 
non-obvious relationship between two classes of mean functionals. Exponential 
tilting is merely a catchy phrase for the operation of applying an exponential 
change of measure to a density or more general measure. In mathematical fi- 
nance and other applications it is known as an Esscher Transform which is a 
key tool for option pricing. We mention that there is much known about expo- 
nential tilting of infinitely divisible random variables and in fact Bondesson (3, 
example 3.2.5) discusses explicitly the case of GGC random variables, albeit 
not in the way we shall describe it. In addition, examining thegamma repre- 
sentation in (jl.lj) one can see a relationship to Lo and Weng (|32l . Proposition 
3.1) (see also Kiichler and Sorensen |28l ) and James |2fj| . Proposition 2.1) for 
results on exponential tilting of Levy processes). However, here our focus is on 
the properties of related mean functionals which leads to genuinely new insights. 

Before we elaborate on this, we describe generically what we mean by expo- 
nential tilting. Suppose that T denotes an arbitrary positive random variable 
with density, say fx- It follows that for each positive c, the random variable cT 
is well-defined and has density 

-Mt/c). 

c 

Exponential tilting refers to the exponential change of measure resulting in a 
random variable, say T c , defined by the density 

_ <r\l/c)f T {t/c) 
W - E[e" cT ] ■ 

Thus from the random variable T one gets a family of random variables (T c ; c > 
0). Obviously the density for each T c does not differ much. However something 
interesting happens when T is a scale mixture of a gamma random variables, 
i.e., T = GgM, for some random positive random variable M independent of 
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Gg. In that case one can show, see (|2ll ). that T c = GgM c where M c is suf- 



ficiently distinct for each value of c. We demonstrate this for the case where 
M = M e (Fx). 

First note that obviously, cMg(Fx) — Mg(F c x), for each c > 0, which in 
itself is not a very interesting transformation. Now setting Tg — GgMg(Fx) 
with density denoted as gg,F x , the corresponding random variable Tg tC resulting 
from exponential tilting has density 

c-'(l/c)j 9!fx ((/c)e^( c ) (3.1) 

and Laplace transform 

V\p-c{l+X)G e M e (F x )] 

E [ e -cG 9 M e (F x )] e • ^' Z ^ 

Now for each c > 0, define the random variable 



(cX+1) 

That is, the cdf of the random variable A c , can be expressed as, 

F Ac (y) = F x (^-^ far0<y<l. 

In the next theorem we will show that Mg(Fx) relates to the family of mean 
functional (Mg(FA c );c > 0), by the tilting operation described above. More- 
over, we will describe the relationship between their densities. 

Theorem 3.1. Suppose that X has distribution Fx and for each c > 0, A c = 
cX/(cX + 1) is a random variable with distribution Fa c - For each 9 > 0, let 
Tg = GgMg(Fx) denote a GGC (8, Fx) random variable having density gg t F x - 
Let Tg >c denote a random variable with density and Laplace transform described 
by i3.1\) and \3.2\l respectively. Then Tg c is a GGC(9, Fa c ) random variable 
and hence representable as GgMg(FA a ). Furthermore, the following relationships 
exists between the densities of the mean functionals Mg(Fx) and Mg(FA a ). 

(i) Suppose that the density of Mg(Fx), say, £#_f x * s known. Then the density 
of Mg{FA a ) is expressible as 

for < y < 1. 

(ii) Conversely, if the density of Mg(F Ac ), £,8F Ac {y)> * s known then the density 
of Mg{Fx) is given by 

2 f n „ ( X \ e-O+rxQ- 



&f x (x) = (1 +x) & Fa . -t—- e 
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Proof. We proceed by first examining the Levy exponent [rpF x (c(l + A)) — 
ipF x ( c )] associated with Tg :C as described in (|3.2p . Notice that 

/>oo 

^ X (c(l + A)) = / Iog(l + c(l + A)a:)Fx(da:) 
Jo 

and ipFxi ) is °f the same form with A = 0. Hence isolating the logarithmic 
terms we can focus on the difference 

log(l + c(l + \)x) - log(l + ex). 

This is equivalent to 

\og(l + -^x\ =logf-J— + _^_(i + A)V 



1 + cx / \ 1 + cx 1 + cx 

showing that Te lC is GGCJj9, This fact can also be deduced from Proposi- 

tion 3.1 in Lo and Weng (|32l). The next step is to identify the density of Mq(Fa c ), 
in terms of the density of Mg(Fx)- Using the fact that Tg = GgMg(Fx), one 
may write the density of Tg in terms of a gamma mixture as 

j.6 — X poo 

9e, Fx (t) = Y^J Q e- t/m ™- e tieF x (™)dm. 

Hence, rearranging terms in (|3.ip . it follows that the density of Tg tC can be 
written as 

.0-1 ^.oo 

Jo 

Now further algebraic manipulation makes this look like a mixture of a gamma(#, 1) 
random variable, as follows, 

cm + 1 



I _ t £E±I 



T(0) 



o 



——————— — C?771 . 

(1 + cm) 



Hence it is evident that Alg(F^ c ) has the same distribution as a random variable 
cM/(cM + 1) where M has density 

e «^x( c )(l + cm )-% Fx ( m ). 
Thus statements (i) and (ii) follow. □ 



3.1. Tilting and Beta Scaling 

This section describes what happens when one applies the exponentially tilting 
operation relative to a mean functional resulting from beta scaling. Recall that 
the tilting operation applied to GgMg(Fx) described in the previous section sets 
up a relationship between Mg(Fx) and M$(Fa c )- Consider the random variable 
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0e*,6{i-<T)Mg a (F x ) = Mg{F XY<r ), Then tilting GgMg(F XYrr ) as in the previous 
section leads to the random variable GgMg{F cXYa / r cXYa+ W) and hence relates 

Pea,e(i-a)M e<T (F x ) = M e (F XY J 
to the Dirichlet mean of order 6, 

M e{FcXY„ / (cXY„ + l)) ■ 

Now letting Fa c y b denote the distribution of A c Y a , one has 

d cXY a 



(cXY a + 1) 



and hence 



M e {F cXYa/(cXY „ +1) ) 1 M fl (*U.O = Pe a ,e{i-<,)Mecr{F Ac ). (3.3) 

In a way this shows that the order of beta scaling and tilting can be inter- 
changed. We now derive a result for the cases of Mi(F XY<r ) = (3 a .i- a M a (F x ) 
and M 1 (Fa cY<j ) = /3 CT) i_ CT M cr (F J 4 c ), related by the tilting operation described 
above. Combining Theorem 12.21 with Theorem 13.11 leads to the following result. 

Proposition 3.1. For each < a < 1, the random variables Mi(F XY<T ) = 
f3 a ,i-<jM a (F x ) and Mi(F AcY „) = P a ,\- a M a {FA e ) satisfy the following; 



(i) The density of Mi(Fa cYcj ) is expressible as 



£,F AcYa {y) 



e aip Fx (c)ya-l 



nF 



XY a 



c(i - y) 



for < y < 1. 

(ii) Conversely, the density of M\ {F XYcr ) is given by 



6 



.(*) 



7T(1 



1 



Proof. For clarity statement [(i)] is obtained by first using Theorem 13. II Which 
gives, 



£,F AcYa (y) 



i 



.(c) 



y 



c(l - y) 



for < y < 1 . It then remains to substitute the form of the density (|2.5[) given 
in Theorem 12.21 Statement [(ii)] proceeds in the same way using (|2.6p . □ 



Note that even if one can calculate <&fa c f° r some fixed value of c, it may not 
be so obvious how to calculate it for another value. The previous results allow 
us to relate their calculation to that of &f x as described next. 
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Proposition 3.2. Set A c = cX/(cX + 1) and define <& Fa (y) = E[log(|y - 
A c \)l{A c ^y)]. Then for < y < 1, 

$F Ac (y) = *f x ( c(1 %) ) - ^ ( c ) + to s(c(l - 2/))- 

Proof. The result can be deduced by using Proposition 13. II in the case of a = 1. 
First notice that smiirFx ( grjzyj )) = sin(7rJ r J 4 c (y)) ■ Now equating the form of 
the density of M\(Fa c ) given by (|1.5p with the expression given in Proposi- 
tion [34] It follows that 

c -»„ (v)_ 6 ^ (C) e -^(^y) 

which yields the result. □ 

Remark 3.1. M^e point out that if G K represents a gamma random variable for 
k, ^ 9, independent of Mg(Fx), it is not necessarily true that G k Mq(Fx) is a 
GGC random variable. For this to be true Mg(Fx) would need to be equivalent 
in distribution to some M K {Fp). In that case, our results above would be applied 
for a GGC(k, Fr) model. We will encounter a variation of such a situation in 
section \4-l\ 



4. First Examples 

In this section we will demonstrate how our results in section [2] and [3] can be 
applied to extend and explain results related to two well known cases of Dirichlet 
mean functionals. 



4-1. An example connected to Diaconis and Kemperman 

Set X — U, denoting a uniform[0, 1] random variable with Fu(u) — u for < 
it < 1. It is known from Diaconis and Kemperman (jlOl ) that the density of 
M^Fu) is 

- sm(Try)y- y (1 - yT (1 ~ v) for < y < 1. (4.1) 

7T 

Note furthermore that T\ = G\Mi{Fjj) is GGC(1, Fjj) and has a rather strange 
Laplace transform, 

We can use this fact combined with the previous results to obtain a new explicit 
expression for the density of what we believe should be an important mean 
functional and corresponding infinitely divisible random variable. 
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Proposition 4.1. Let W = G\/G'i be the ratio of two independent exponential 
(1) random variables having density F\y{dx)/dx = (1 + x)~ 2 for x > 0. Now 
let Ti denote a GGC(1, random variable, with^p w (\) = -^j-log(A). Then 

T\ = GiMi(Fw), where M\(Fw) has density, 

1 / ttx \ x 

6v [x) = - sm — — x its /or ac > 0. 

7T \ 1 + X J 

Proof. First note that it is straightforward to show that 

E[e~ Tl ] = E[(l + MxiFw))- 1 ] = e" 1 . 

This fact also establishes the existence of T\ . Now we see that 

W G t d 



W + l (Gi + G'i) 



U. 



Hence Theorem l3.1l shows that M 1 (Fjj) arises from tilting the density of GiMi(Fw) 
The density is obtained by applying statement (ii) of Theorem 13.11 to (|4.ip . Or 
by statement (ii) of Proposition [37X1 □ 



We now apply Theorem 12.21 and Proposition 12.11 to give a description of 
the finite-dimensional distribution of the subordinators associated with the two 
random variables above. Here U and W = G\/G\ are as described previously. 

Proposition 4.2. Let U denote a uniform^}, 1] random variable and let W = 
G\/G'i denote a ratio of independent exponential(l) random variables. 

(i) Suppose that (£i(i);0 < t < 1) is a GGC{\, Fjj) subordinator, then for 
(Ci, . . . , Ck) a disjoint partition of (0, 1), the finite- dimensional distribu- 
tion has a joint density as in \2. 7\ ), with, 

g^Fvixi) = l\-^—^a l {l-y))y^ 1 'y^{l^y)-' ydl ^ ) dy, 
Jo k 

for i — 1, . . . , k. 

(ii) That is Ci(Ci) = G\Mi{FxjYv.) and are independent for i = 1, ...,fe. 
Furthermore, the density of Mi(Fjjy„.) = fioi.i-aiMa^Fjj) is 

— sin(^(l - v))v°*W-\l - y)-°&-ti 

IT 

for < y < 1. 

(Hi) //(Ci(^);0 < t < 1) is a GGC(1, Fyy) subordinator then the finite- dimensional 
distribution (Ci(Ci), . . . , Ci(C&)) * s described now with 

9c % .F w {xi) = / e~^-sm( ) w^~ 3 dw. 

,/n 7T V 1 + W J 
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(iv) That is £i(Ci) = G\Mi{Fwy„.) and are independent for i = 1, ...,k. 
Furthermore, the density of M\{FwY a .) = lia i .\-a i M ai [Fy/') is 



— sin I — ] .i' 1 

7T V 1 + X , 



-1 



for x > 0. 

Proof. This now follows from Theorem 12. 2\ Proposition 12.11 and (|4.ip . Specifi- 
cally, note that C(Fu) = (0, oo), hence for any < a < 1, 

sin(7ri 7 j/Y'„ (w)) = sin(7rcr(l — «)) 

for < it < 1 and otherwise. Furthermore from (|4.1[) . or by direct argument, 
it is easy to see that, 



^(y) = -log (y-Hl-y)- (1 - v) ) 



1. 



This fact also is evident from Diaconis and Kemperman fiol ). It follows that 
Mi {Fuy„ ) has density 

— sin(7ra(l - y))^ 1 ^" 1 ^ - y) _<T(1_a) for < y < 1. 

The density for Mi (Fwy„ ) is obtained in a similar fashion by Proposition 13.11 

□ 

Remark 4.1. Setting 

d cGi 



cGi + G'i ' 

one can easily obtain the density of the random variable Mi (F^ c ) for each c > 
by using statement (ii) of Theorem \3.1\ Note also that one can deduce from the 
density of Mi(Fw) that $>f w {x) = + ^)]log(^)- Hence in this case an 

application of Proposition \3.°2\ shows that. 



c(l-y)+y 



y 



~(i - y) 



clog(c) + lo g[c (l- y )]. 



c- 1 



We note that otherwise it is not easy to calculate ^a c , in this case, by direct 
arguments. 



4-2. Reconciling some results of Cifarelli and Mellili 

To further illustrate our point we show how to reconcile two apparently unre- 
lated results given in Cifarelli and Melilli (d). Let /?i/2.i/2 denote a beta(l/2,l/2) 
random variable with cdf 
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often referred to as the arcsine law. Cifarelli and Melilli (8, p. 1394-1395) show 

that for all d > 0, Mg(Fp 1/21/2 ) = (3 e+ 1/2,6+1/2- Now define the probability 
density 

Qi/2{dx)/dx = — a: _1 / 2 (l + x) , 

Cifarelli and Mellili (d, p.1394-1395) also show that for 6 > 1, M e (g 1/2 ) has the 
density, 

Hjort and Ongaro flih recently extend this result for all 9 > 0. 

Here, however, we note that if X has distribution g 1 / 2 then X = Gi/ 2 /G^ 2 
where G1/2 and are independent and identically distributed gamma ran- 
dom variables. That is Mg(gi/ 2 ) = Mg(Fa 1/2 /G' 1 2 )< Now using the known fact 
that 

d G 1/2 /G' 1/2 
Pi/2,1/2 = 



Gi/ 2 /G' l/2 + 1 

we see that /?i/2,i/2 is a special case of ^! in Proposition ^. 21 That is, Mg(Fp 1/2 ) 
results from exponentially tilting the density of GgMg(F Gl/2 / G ^^) 1 with c = 1. 

Hence one could have obtained the density of Mg(FQ 1/3 /G> ) by using the re- 
sult of Cifarelli and Melili that Mg(Fp 1/2 ) = /?e+i/2,e+i/2j an d simply applying 
statement (ii) of Theorem 13.11 

Remark 4.2. Note that the following identity for the normalizing constant, 

r(20 + i) _ 4 _ e= r(e + i) 



r(0 + i/2)r(0 + i/2) t(6 + i/2)r(i/2)' 

is obtained from the gamma duplication formula, where in particular 

T{29 + 1) = W/2) T{e + 1)T{9 + 1/2) - 

Hereafter we shall use the notation Xi/ 2 = Gi/ 2 /G' 1 ^ 2 . Moreover, from the 
density in (|4.2|) . it follows that, using the notation Xing, 

Y — M (TP \± G S+l/2 
^1/2 

In particular Xi/2,0 = -^1/2- Furthermore, 

d d Gg 1 2 
GgXi/ 2t g = GgMg(Fx 1/2 ) — Gg—— 

(-a/2 

is GGC(9, Fx 1/2 ), with Laplace transform 

E[e- AG9M9(F *v 2 )] = (1 + ^)" 2e . (4.3) 
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Now, for < p < 1, setting c = p 2 /(l — p) one can extend the result of Cifarelli 
and Melilli for mean functionals denned by the random variable, 

_ P 2 Gi /2 _ cXl/2 

1/2 ' p ~ {p 2 G 1/2 + q*G> 1/2 } ~ cX 1/2 + l 

having density 

qp y-^jl - y)~ 1/2 

2 — I — 27i r for < y < 1, 

7T q z y+p z (l-y) 

as follows, first the density of the random variable 

cM e (F Xl/2 ) = M e (F cXl/2 ) 

is simply, 

CeFx {x ) = F ( g + X ) ci/V-^Cc + x)-^. 

w«xi /a v ; r(6> + |)r(|) v ; 

The density of Mg(Fs 1/2 p ) is then obtained by applying statement (ii) of The- 
orem [3J] to the density above, to get, 

r(g + l) y e -VHl-y) e - 1/2 

6 ^ W=W ^7R+?F" (4 ' 4) 

However this density equates with the density of a two parameter (1/2,0) 
Poisson-Dirichlet random probability measure evaluated at a sect C say, -Pi/2,e(C)> 
such that E[P 1 / 2 ,e(C)] = p, given in Carlton (0, Remark 3.1) (see also |23|) 

and (34)). Hence we see that M g {F Bl/2p ) = P 1/2 ,e (C) arises from exponentially 
tilting the density of GgMg(Fx 1/2 ). 



4.2.1. Results for 9 > -1/2 
Notice that the random variable 

v d G e+1 / 2 

A l/2,8 - ~~F< 

^1/2 

is well defined for all 9 > —1/2, and still has density (|4.2p for the negative 
range of 9. However obviously in the range — 1/2 < 9 < it is not a Dirichlet 
mean functional of order 9, and as such is not covered by the result in (§). The 
next result shows that for all 9 > —1/2 the random variable Xi/2,8 is a mean 
functional of order 1 + 9. 

Proposition 4.3. For each 9 > —1/2, the random variable Xi/2.6 satisfies the 
distributional equality. 

Xi/2,8 = ^9+1/2,1/2^1/2,9+1/2 = Ml+e(Fx 1/2 Y e+1/2 ) 

1+0 
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As a consequence the random variable, G\+gXi/2,e — Ge+i/2^i/2,e+i/2 *s GGC(9+ 
1) ^y 1/2 y 9+1/2 ) =GGC(9 + 1/2, Fx). 

1 + 

Proof. To see this note the simple identity 

Gg+l/2 = ^9+1/2,1/2^9+1- 

The first equality in the Proposition is then obtained by dividing both sides by 
an independent Gi/ 2 random variable. The second equality is obtained by using 

TheoremEHJ That is using the fact that e+ i/ 2 , 1/2^1/2,0+1/2 = A?+i/2.i/2^0+i/2(-Fx 1/2 )■ 

□ 

Thus we see that for 6 > 0, 

Xi/2.6 = Mg(F Xl/2 ) = Mi + g(F Xl/2 Y e+1/2 )• 

Remark 4.3. The random variable = Gi^/G^^ is a special case of an 

important random variable studied by Lamperti K291 ). The general case is de- 
fined for < a < \, X a = S a /S' a> where S a and S' a are iid positive stable 
random variables of index a. The GGC(9, Fx 1/2 ) random defined by its Laplace 
transform {4-.Sfy is an example of the generalized positive Linnik random variable 
which can be represented as (Gg/a) 1 ^ 01 S a . Where in this case a = 1/2. A thor- 



ough investigation of this class, and its various implications, initiated in \2& ). 
will be discussed in a forthcoming manuscript YE$i ) (see also \EK)) 



5. Obtaining the finite dimensional distribution of a subordinator of 
BFRY 

The previous examples revisited some cases that have existed in the literature 
for some time. Our final example shows how one can apply the results in section 
2 to obtain new results for a subordinator recently studied by Bertoin, Fujita, 
Roynette and Yor (0). We first supply some background. For < a < 1, let S a 
denote a positive a-stable random variable specified by its Laplace transform 

E[e" AS °] = e" A< \ 

In addition define 

z a = (f ) 

where S' a is independent of S a and has the same distribution. The density of this 
random variable was obtained by Lamperti(29)(see also Chaumont and Yor ((0), 
exercise 4.2.1) and has the remarkably simple form, 

sin(7ra) 1 

ira y z + 2ycos{na) + 1 
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In order to avoid confusion we will now denote relevant random variables 
appearing originally in (Q) as A a and G a , as S a and G a respectively. From (0), 
let (£ Q (t); t > 0) denote a subordinator such that 



E [ e -AE a (t)] _ ((A + 1)«-A a ) 

= exp(-t(l - a)E[log(l + 1/G a )]) 

where from (Qj, Theorems 1.1 and 1.3), G Q denotes a random variable such 
that 



and has density on (0,1) given by 

asin(7ra) u a_1 (l — u) c 



(1 - a)7r u 2q _ 2(1 - u )a cos(7ra) + (1 - u) 



2n 



Hence it follows that the random variable 1/G a takes its values on (1, oo) with 
probability one and has cdf satisfying, 

1 - F 1/Ga (x) = F Ga (l/x) = F Zl _ a ((x - 1)- Q ), 

where using properties of the cdf Zi- a , (see tflTh and ( (|23 ). Proposition 2.1)), 
it follows that for x > 1 



^ ( 1 ) = i _ 1 co t-i " «)) + (» - 1 )" 



'(x— 1) 7r(l — a) \ sin(7r(l — a)) 

It follows as noted by (0) that (£ Q (i); f > 0) is a GGC(1 — a, F 1 / Ga ) subordi- 
nator. Where from (|2), Theorem 1.1) the GGC(1 — a, F 1 / Ga ) random variable 

E Q = S a (l) satisfies, 

where U denotes a uniform[0, 1] random variable and for clarity Gi_ Q is a 
gamma(l — a, 1) random variable. It is evident, as investigated in l)l7l). that 

Afi_ a (Fi /G J = -3^- = 

Pa,l 

So this constitutes a rare known example where a Dirichlet mean of order < 

6 < 1 has a simple recognizable form. Note also that Mi_ a (F Ga ) = U for all 
< a < 1. The above points may also be found in the survey paper of (|24T ). 

Note that when a = 1/2, G1/2 = /3i/2,i/2> which then relates to the results 
in section 4.2. In addition, it is known that for each fixed t, 



d G t / 2 



/3l/2,(l+t)/2 
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This result may be found in (|25l ) where the case of a = 1/2 also belongs to 
another class of subordinators indexed by a. 

Now using the fact discussed in this section we will show how to use the 
results in section 2 to explicitly describe the finite dimensional distribution of 
the subordinator T, a (t) over the range < t < 1/(1 — a), hence by infinite 
divisibility for all t. Additionally the analysis will also yield expressions for 
mean functionals based on F 1 / Ga . We first show how to calculate 

n a (x) = $ Fi/Gq (x) = Epogflt - l/G a |)I(t + 1/G Q )]. (5.1) 

Using simple beta gamma algebra one has 

^ d_ Gl—a d_ *i Pi — a. a 
— ~7, — "I - 



(3 a ,i C/ 1 /- 

Hence applying Theorem 2.1, with 0—1, and a = 1 — a, it follows that S Q is 
also GGC(1, Fyi-a/Ga) and 

Pi— a a d 01— a a d . , , „ , /r o \ 

= = Mi(F Yl _ a/Ga ) (5-2) 

Furthermore, using (|2.4| . one has 

• / v r \\ / sin(7r(l - a)Fa (1/x)), if x > 1, 

sm 7riv, /G (#)) = > • 1 rt w -r r. ^1 (5-o) 

v n-a/ibrav. /; ^ sm(7r(l - a)), if < x < 1. v y 

where again using the properties of Fz x _ , as deduced from f (|22l ). Proposition 
2-1, [(iii)D, 

^(1 - a)F Ga (1/x)) = Sln( " (1 " a)) m (5.4) 

[{x - 1) 2Q - 2(x - l) a cos(Tra) + 1] 

for x > 1. These points lead to the following description of (|5.1[) . 

Proposition 5.1. For < a < 1, consider lZ a (x) as defined in h5.1\) . Then, 

n a (x) = I w [1 r ^ r f-vr" 2( ^A )aco8(7ra)+ij)] ' ?! J \ 1 ^ 

i T^^S^/t 1 ~ l 1 ~ x > D' «/0<x<l. 
Proof. By standard calculations the density of B a — Pi- a , a / Pa,i is given by 

f Ba (x) = Sin( " (1 ~ Q)) x— Ml - (1 - < 1)] 

However we see from (15.21) that £? Q = Mi(i ; y 1 _ /6 0[ )' Hence Theorem 2.2 applies 
and the density of B a can be written as 

f B Jx) = ?— s i TL (irF Yl _ a/Ga (x))e-( 1 - a ^'> 

Now equating the two forms of the density of B a and using (|5.3p and (I5.4[) , one 
then obtains the expression for lZ a □ 
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This leads to our next result, 
Proposition 5.2. For < a < 1 the density of Mi(F Ya /G a ) is given by, 

fa,a{x) = Sm(7T(T)[l - (1 - x) a }~ 

7T 

for < x < 1, and for x > 1 

2^ I 

x i— 



/a,a(ic) = sin(^F Go (l/a;))[(a; - l) iQ - 2{x - l) a cosfra) + 1] 



7T 



Proof. From Theorem 2.2 we have the general form of the density of Mi(Fy cr /& a 
is given by 



sm(7rF Ya/Ga (x))e-^^ 



7T 

The result is then concluded by applying Proposition 5.1 and (|5.3[) and (|5.4[) . □ 

We now can close with the finite dimensional distribution of the subordinator, 
which follows from Proposition 2.1 and Proposition 5.2. 

Theorem 5.1. Consider the GGC{1 — a, i*i/G a ) subordinator (E Q (i),i < 1/(1 — 
a)). XTien /or a disjoint partition of the interval (0, 1/(1 — a)], (Ci, . . . , Cfe), £/ie 
/imie dimensional distribution o/(£ Q (Ci), . . . , £ a (Cfe)) is swc/i iftai eac/i £ a (C;) 
is independent with density 

e -*/™ w -lf aai ( w )d w 

where a i = (l — a)\Ci\ and f atai (w) represents the density of Mi{Fy a . /g„) given 
in Proposition 5.2. That is, 

^ a {Ci) = GiMi(F Yr ./G a )- 
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